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ABSTRACT 


We derive the source functions and the energy fluxes for wave generation in magnetic flux tubes embed- 
ded in an otherwise magnetic field-free, turbulent, and compressible fluid. Specific results for the gen- 
eration of longitudinal tube waves are presented. 
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1. INTRODUCTION 

Wave generation by turbulent motions in the outer convection zones of stars has long been thought 
to be central to the heating of stellar chromospheres and coronae (e.g., Biermann 1946 and 
Schwarzschild 1948). The early suggestions were followed by number of detailed studies in which the 
generation of acoustic waves (Lighthill 1952, Proudman 1952, Stein 1967, Renzini et al. 1977 and Bohn 

1980, 1984) as well as MHD waves (Kulsrud 1955, Osterbrock 1961, Parker 1964, Kuperus 1965, Stein 

1981, Ulmschneider and Stein 1982, Musielak and Rosner 1987a, b) was considered. These latter calcu- 
lations are all based on the assumption of a uniform and weak background magnetic field, an assumption 
which is contradicted by solar and stellar observational evidence for inhomogeneous and locally strong 
magnetic fields (cf., Harvey 1977, Stenflo 1978, Robinson et al. 1980); thus, at least the solar magnetic 
field has instead a "flux tube" structure, and flux tube waves carrying the wave energy away from the 
convection zone may well be responsible for heating at least some portion of the outer atmospheric 
layers (cf.. Spruit and Roberts 1983). In addition, it has been shown that the acoustic and MHD energy 
fluxes generated in stellar convection zones - based on calculations which assume homogeneous mag- 
netic fields - are insufficient to explain the UV and soft X-ray fluxes observed by the 1UE and Einstein 
Observatories (Linsky 1981, Vaiana et al. 1981, Ulmschneider and Bohn 1981, Rosner et al. 1985, 
Musielak and Rosner 1987b). 

In this paper, we derive the source function and energy flux for flux tube waves. We consider 
magnetic flux tubes embedded in a magnetic field-free, turbulent and compressible medium, and assume 
that the tubes are thin and oriented vertically, the latter assumption allows us to separate the generation 
of compressional tube waves from incompressional waves. In the present paper (Paper I), we concen- 
trate on the generation of longitudinal tube waves, discuss the wave propagator and the relevant critical 
frequencies, and finally discuss the dependence of energy fluxes on the parameters which enter into the 
calculations. The generation of transverse tube waves and the application of the results to late-type stars 
will be treated in following papers. 

The plan of our paper is as follows: The MHD equations and the basic formulation are presented 
in Section 2; the inhomogeneous wave equation and its solutions are given in Section 3; the energy 
fluxes for longitudinal tube waves are described in Section 4; and the model parameter dependence of 
energy fluxes and their discussion are to be found in Section 5. A summary of our new results and our 
conclusions are given in Section 6. Two appendices contain mathematical details, which amplify discus- 
sions in the main text 
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2. MHD EQUATIONS AND THE BASIC FORMULATION 

In this section, we discuss the basic equations of motion, and develop the formalism for calculating 
the generation rate for magnetic tube waves. In order to simplify the problem to the essentials, we shall 
assume that the fluid is locally isothermal, that the gas pressure is a scalar, and that displacement 
currents and electrostatic forces may be neglected; furthermore, it is straightforward to show that (as 
long as shock formation does not occur) dissipation by molecular viscosity and Ohmic diffusion is negli- 
gible for the problem at hand. In the following, we present the linearized magnetohydrodynamic (MHD) 
wave equations, the basic assumptions for flux tubes, and finally the set of equations used to calculate 
the rate of tube wave generation. 


(a) The Linearized MHD Equations 

Our assumptions lead to the ideal MHD equations, which may be written in the following linear- 
ized form: 


|-P + V-(P.tf> = N t , 


( 2 . 1 ) 
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( 2 . 2 ) 


and 


J-* 

dt 


[(Vx^xjf, ] = 


(2.3) 


j;# - Vx(*x#,) = AV. 


(2.4) 


where p 0 ,p 0 ,B 0 refer to the unperturbed atmosphere, p, p, U, 2? are the perturbations of density, pres- 
sure, velocity and magnetic field, respectively, and V, ( a [*yR77p] 1/2 ) is the sound velocity; the magnetic 
field is assumed to be potential, and the equation of motion is simplified by the assumption that the 
background atmosphere is in static equilibrium. 

All terms linear in the perturbations in equations (2.1) - (2.4) are written on the LHS; these terms 
define the wave propagation operators for MHD waves. However, terms quadratic in the perturbations 
are collected on the RHS, and are treated as known quantities described by a given flow (Lighthill 1952, 
Stein 1967). The latter terms determine the source function responsible for the MHD wave generation 
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(Musielak and Rosner 1987a), and are defined as follows: 


Ni = " V (pit ) , (2.5) 

N 2 = - i? • V p + V 2 it • V p , (2.6) 

#3 = -Sr 4: * -■£-<*• V)J? + _L_ (Vx^xB*. (2.7) 

P« of p 0 4jrp 0 

= Vx^xfi*). (2.8) 


(6) 77ie 77u'n F/ux 7u£e Approximation 

In order to calculate the energy fluxes for flux tube waves, we assume that a vertically-oriented 
magnetic flux tube is embedded in a stratified and non-magnetized medium, and that all unperturbed and 
perturbed quantities depend on z and t alone; the z-axis is identified with the vertical direction (i.e., 
g = - gz ) which, in our approach, is just the tube axis. We then obtain 


p = p(zj) + 0(e), (2.9) 

p = p(zj) + 0(e) , (2.10) 

it • u,(zj) 2 + 0(e) , (2.11) 

$ = B,(zj)z + 0(e), (2.12) 


where e is defined as the ratio of the tube radius R, to the tube length L,; for e < 1, the tube can be 
treated in the thin flux tube approximation, an approximation which allows us to consider all unper- 
turbed and perturbed quantities to zeroth-order (see, for example, Roberts and Webb 1978). In this 
approach, the magnetic field within the flux tube is essentially axial, and is described by - B 0 (z)z; 
note that the solenoidal condition does not restrict B 0 (z), but does however allow one to calculate the 
horizontal components of the magnetic field once the vertical component is known. 



• + I 


-5 - 

We further assume that the cross-section of the tube is always circular, and that the thin flux tube 
is in temperature equilibrium with the surroundings, so that both density gradients inside [p„ ( 2 )] and out- 
side [p,(z>] the tube are described by the same density scale height; hence, the vertical density variation 
is given by 


Po*b) 


Poo exp 



(2-13) 


where ( a V, 2 iyg) is the density scale height, and is identical to the pressure scale height for an isoth- 
ermal atmosphere. The vertical magnetic field variation is given by 


B 0 b) - Boo exp 



where the magnetic scale height H b can be defined as 



Po 

P« - Po 


= 2 H 


p > 


(2.14) 


(2-15) 


and where V A 2 ( 3 B?l4np a ) is the square of the Alfven velocity. 

In addition, if the total (gas plus magnetic) pressure is constant across the tube, then we must have 
(see Roberts and Webb 1978) 

P.b) = Pob) + ~~ B?b ) , 

and 

B?b) _ Pob) = p,b) 

Bo\ 0) " PoiO) pA 0) ’ 

which allows us to easily calculate the structure of the tube when the boundary values for the pressure 
and magnetic field are given. Finally, we assume that the tube is untwisted, so that we ignore processes 
— such as the generation and propagation of torsional tube waves - as well as instabilities such as the 
lrinic instability. Such transverse wave motions will be considered in a subsequent paper. 


(2.16) 


(2.17) 
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(c) The Flux Tube Equations 

Combining the continuity and induction equations, and using equations (2.9) - (2.12) and (2.13) - 
(2.15), the set of MHD equations (2.1) - (2.4) can be re-written as 


ar* 


-£*> 


+ W 2 U, = It! , 


(2.18) 


where n x contains the nonlinear terms, and is given by nj = ~{\lp 0 )N x - (1 IB 0 )N 4 ,\ 

^ Pi ~ V, 2 -jj- pi + W g u, = n 2 , 


(2.19) 


where n 2 = (1/p olW* 


d 

— u, + W ipi + pig = ji 3 ; (2.20) 

where n 3 = N 3t ; N 2 , and N 4 , are the z-components of ^ 3 and (equations 2.7 and 2.8), respectively. In 
addition, p x = p/p ot p\ = pip a , and B x = B,IB 0 are the new dimensionless perturbations; the operators W x , 
W 2 and W t are defined by 




3_ 

dz 
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w, * 4- - 


3_ 

dz 


1 

2H 0 


W. 


( 1 ~ Dg 


( 221 ) 


The set of equations presented above (plus the horizontal pressure balance equation; see next sub- 
section) fully describes vertical motions inside the flux tube; The non-linear terms n u n 2 and n 3 can be 
treated as the source of these motions, and can in principle be determined once the driving turbulent 
flow is specified. Such an approach has been followed by Stein (1968) for the purely acoustic case, and 
by Musielak and Rosner (1987a) for the case of homogeneous magnetic fields. As discussed below, we 
assume in our present (strong magnetic flux tube field) case that these internal sources are suppressed by 
the internal magnetic field, so that the wave motions inside the tube are driven from the outside alone. 
If there are no turbulent motions within the tube, the equations describe the propagation of longitudinal 
waves along the tube (see Roberts and Webb 1978 for comparison), and can be easily solved since all 
coefficients in these equations are constant To close this set of equations, we need a relationship 
between the total pressure (gas and magnetic) inside the tube and the external pressure. 



(d) The Horizontal Pressure Balance 

The relation between the variable fluid pressure applied to die undisturbed tube boundary by the 
external turbulence (or any other pressure variations) and the internal fluid pressure can be written as 
(Parker 1979) 


p, (*J) + -^5/W) = P.(*J), (2.22) 

where p, and B, are the gas pressure and the magnetic field inside the tube, respectively. The pressure 
P t outside the tube is defined by 


P* - P.(z) + + Pu(*J), (2-23) 

where p, is the gas pressure of the external fluid in the absence of any turbulence; Purb(zj) is the pres- 
sure of the external turbulence and p u (zj) is the disturbance of the external pressure due to the back- 
reaction of the moving tube. The last term disappears if the tube is at rest with respect to its surround- 
ings; this term will be neglected in our calculations. Linearizing equation (2.22) and using equation 
(2.16), one finds 


Px + VZB X 



l 

Pturb • 
Pff 


(2.24) 


where the variable fluid pressure pa^ caused by the external turbulence can be defined as 


Pturb 



(2.25) 


and where u, is the turbulent velocity. 

Using the equipaitition energy to calculate Bf, and with the help of equation (2.23), one obtains 

Px + Va Bx — » (2.26) 


where 
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and where u a and are the x- and y-components of the turbulent velocity, respectively. 

The set of tube equations (2.18) through (221), together with equation (226), fully describe the 
generation (due to the turbulent motions outside and possibly inside of the tube) and the propagation of 
longitudinal tube waves; this set of equations can be written as a single inhomogeneous wave equation, 
and can be solved by Fourier-transforming this equation. This problem will be addressed in the next 
section. 
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3. THE INHOMOGENEOUS WAVE EQUATION 


In this section, we derive and solve the inhomogeneous wave equation for longitudinal tube waves. 
We also calculate die source function, and discuss the cutoff frequencies for longitudinal tube waves. 

(a) The Wave Propagator 

We eliminate the density, magnetic field and velocity perturbations from equations (2.18) - (2.21) 
and (2.26), and obtain the inhomogeneous wave equation for the pressure perturbations in the form 


JL _ v *jL 

dt 2 ‘ dz 2 


+ 2 V, ( oo* - <a bl ) ~ + 


P, a>?t 


Pi = S p (zj), 


(3.1) 


where the characteristic velocity for longitudinal tube waves (Defouw 1976) is given by 


and 


V, 2 = 


Vt V , 2 

V 2 + V 2 ' 



1 + 



(3.2) 


(3.3) 


we also define the three critical frequencies o) a , co*. and co^, 




V, 


2 H e 


CD* 


V, 

2 H b ' 


co« 


V, g 

— JL /v - 


V, 


Va V. 


(Y - D“.= -r~ • 

Va 


(3.4) 


where CD*, is the Bnmt-Vaisala frequency. Note that for intense magnetic flux tubes (V, < V A ), c o ct 
becomes the magnetic Brunt-Vaisala frequency, however, for V, •> V A , the characteristic tube velocity 
V, - V A and co a = cot,. 

In equation (3.1), all the non-linear terms are collected on the RHS, where they become the source 
function for the longitudinal wave generation. We postpone the calculations of these source terms to the 
next subsection, and provide here only the definition of the source function 



g *2 



(V 2 n x + n 2 ) + 


Vs_ 

Va 





S p (*4) = 
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+ ( V , 2 w 2 + w t ) 



f * 

a 

n 3-g 

dr 



(3.5) 


To eliminate the first-order space derivative from the inhomogeneous wave equation (equation 3.1), 
we make the following transformation 


Pi - P 2 



(3.6) 


and obtain the inhomogeneous wave equation in the form 




+ n 2 


Pi = S,(zj) , 


where 


S,(z,t) = 



r 

S p (z4) 


and 


£1? — ( (Ogt - Oj, ) 2 + p, (o 2 t ; 


(3.7) 


(3.8) 


(3.9) 


with help of equation (2.15), one obtains instead 


£2, 2 = j co£ + P« co* . (3.10) 

Fmally, after some algebraic manipulations, equation (3.10) can be written in the form given for 
the first time by Defouw (1976) 


Zl _L _ _L YL 1 - 1 

16 2y V} -f . 


(3.11) 


The LHS of equation (3.7) is the propagation operator for the tube waves; this relation allows for 
cross-sectional variations. In this case, the gas pressure is the principal restoring force. The waves 
which result are longitudinal tube waves, which may be viewed as acoustic waves propagating along the 
tube, but modified by the tube geometry. As a result of this modification, the characteristic phase velo- 
city equation (3.2) and the critical frequency for die vertical propagation (equations 3.10 and 3.11) are 
different from those for acoustic waves. By studying the dispersion relation (which is global in the 
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approach considered here), it can be shown that ft, is the cutoff frequency for longitudinal tube waves, 
so that the waves cannot propagate if their wave frequency is either lower than or equal to this cutoff. 
As shown by equation (3.10), the cutoff frequency for longitudinal waves is not as simple as for acous- 
tic waves, and depends on both density and magnetic field scale heights, as well as on the Brunt-Vaisala 
frequency modified by the tube geometry, note, however, that the value of the critical frequency ft, is 
always comparable to the acoustic cutoff frequency (see equation 3.11). The tube cutoff frequency 
reaches a maximum when the tube structure is entirely dominated by the magnetic pressure (e.g., when 
the gas pressure inside the tube is negligible); in this limit, however, the longitudinal tube waves cannot 
propagate (this point will be further discussed in section 5b). 


(b) The Source Function 

To calculate the source function defined by equation (3-5), we make one additional assumption, 
namely that there are no turbulent motions inside the thin flux tube; this Ansatz simplifies our calcula- 
tions substantially since = n 2 = n 3 = 0. Under this assumption, it is only the external turbulence 
which is responsible for the wave generation; this can be justified if we recall that strong magnetic fields 
(such as those characterizing the thin flux tube we are considering) will tend to suppress turbulent 
motions within the tube. Equation (3.8) can then be written in the following form 


S,(zj) = 


1 


2(p 0 B e ) 


in 


V, 

2 

' a 2 

_ 0 o 

± + J_‘ 

- m? 



3/ 2 

L 5 

dz H b 



P ,(M m + M„ ), (3.12) 


where M a s u£ and s u^. 

We now expand the source function S,(z,t) in multipoles (Unno 1964), and obtain 


S'izj) = 


2 (fl.) 


in 


Pi 


in 


P ' 2 dt 2 ~ y Pf2 V ‘ dz " 


P« ( ). (3.13) 


where p, a V,IV A and P, a V t IV A . 

The source function does not depend on the second derivative of or M n (no quadrupole emis- 
sion!), but does depend, however, on the first derivative with respect to z (dipole emission) and the 
monopole source term (©«), as well as on the time-dependent term. This latter term becomes negligible 
for stationary or quasi-stationary turbulence. It is interesting that the source function does not depend on 
the cutoff frequency for longitudinal tube waves, and instead depends on both tube critical frequencies 
to* and (o a . The latter fact has important consequences for wave generation (see section 5c), and distin- 
guishes tube wave generation from acoustic wave generation (Stein 1967); for purely acoustic waves, 
both the emission source and the propagation operator depend on the same critical frequencies (the 



- 12 - 


acoustic cutoff and the Bnmt-Vaisala frequency). 

(c) The Solution 

Having obtained the source function, we Fourier-transform equations (3.7) and (3.12) in one 
dimension, using 


[>2(2,0, S,(z.O] = j j dk’d to' \p 2 (k',ai'). S t (i fc'.to')] exp[ HuTt - k'z) ] . (3.14) 

where 

[p 2 (k ,< b ). S t (k ,co)] = jjdz'dt' [p 2 (z V 0 . S, (z 7 0 ] exp[ —i (cor ' - kz')], (3.15) 


and obtain the solution for the pressure perturbations emitted by the turbulent motions. 


p 2 (k,a>) = - 


S,(k, to) 


to 2 - Q, 2 - k 2 V, 2 


(3-16) 


where S, (k ,co) is the Fourier transform of the source function given by equation (3.15). To calculate the 
explicit form of the source function S,(i,< a), we integrate equation (3.15) by parts, and obtain 


tl/2 


S,(k.i o) = P. 


\B 0 


jjdz'dt' 


P, 2 CD 2 + ~ifc P/co^ V, + CO, 


x(M a ' + M„’)exp[-aw'~kz')]. (3.17) 

The solution of the inhomogeneous wave equation given above will be used to calculate the energy 
fluxes emitted as a result of the interaction between the turbulent motions and the flux tube. Note that 
for S,(k,(o) = 0, the solution given by equation (3.16) leads to a dispersion relation which describes the 
propagation of longitudinal tube waves, and the critical frequency £2, becomes the cutoff frequency. 
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4. ENERGY FLUXES FOR LONGITUDINAL TUBE WAVES 

In this section, we derive the final expression for the energy fluxes emitted as longitudinal tube 
waves. We begin with the definition of energy flux for the flux tube geometry and then present and dis- 
cuss the final results. 


(a) The Mean Energy Flux 

The energy flux can be calculated by using the energy conservation principle, which gives 

<F{zjt)> = <pu, A, (|a, 2 + W)>, (4.1) 


where A, is the cross-sectional area of the tube and W is the specific enthalpy; note that the mean energy 
flux defined above is multiplied by the cross-section of the tube, and therefore has dimension [ ergs r -1 ]- 
It should also be noted that within the thin flux tube approximation, the magnetic terms in equation (4.1) 
cancel one another, and that the energy flux does not depend on the magnetic energy; this reflects the 
fact that the gas pressure is the principal restoring force. 

Expanding in a perturbation series, and considering second-order quantities only, one obtains 


< F(zj) > = ~A 0 
2 Po 


<P«*> + ^ 0 <Aa»> + A 0 <pu,> , 


(4.2) 


where A is the perturbation of the tube cross-section. The first term in equation (4.2) can be neglected 
as there is no net flow through the boundary of the tube. One may estimate both remaining terms in 
equation (4.2) using the solutions given by Herbold et al. (1985), 


_5 P.M _5 Xl_ _5_ b2 

2 p«,A e " 2 yV? 2y P ' 


(4.3) 


Finally, by combining equations (4.2) and (4.3), one obtains 


<F(zj) > = A„ (1 + 5) <put> , 


(4.4) 


where 8 = 5P//2y. 
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Using the relation between the pressure perturbations p and p 2 (equation 3.6) and replacing u, by 
Uj*, the energy flux can be expressed in terms of p 2 and written in the form 


<F(zj) > = (p 0 B e ) xn A 0 (1 + 5) <p 2 K»’> . 


(4.5) 


where the velocity perturbation u,* is calculated from equations (2.18) - (2.21) and (2.26), and is given 
by 


where 


m 



1 + toi 




-i 



Pi . 


X a 


4 - y 1 

4y h 9 ■ 


(4-6) 


(4.7) 


Now, we replace «,* in equation (4.5) by (4.6) and obtain 


<F(zj)> = -0,(1 + 8) < Pi 





1 + (Ob. 


± 

dt 


P 2 > 


(4.8) 


where O, ( a A g B 0 ) is the magnetic flux. 

To evaluate the mean energy flux given by the above equation, we express Pifz, t) and pl(zj) in 
terms of its Fourier transform (see equation 3.14), and obtain 


<F(zj)> = 0,(1 + S)jjjjdk'dk"d<o'd(o" a>" x 

X Pi(k',<0') Pi(k",a") exp[j(0)' - <o") t - i(k' - k") z ] . (4.9) 


We take the time average of the mean energy flux by performing die integration over time T 0 , and then 
obtain 


2 * 0 , 

To 





L,(.z&) a <—<F(z^)>>r # 


d + 5) 


Pi /„(z,£D) . 


(4.10) 
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where 


i„{k .CW) = JJdfc' dk" ( k" - a ) expH( *' - *"> 2 1 • (4-1 !) 


P£ = 


©* 


* - - “. 2 


(4.12) 


w 2 - ©£ ~ V , 2 

and where refers to a temporal average over a time scale T 0 . Note that L,(*,©,z) does not have the 
standard energy flux dimension, but rather has dimensions [ergs s' 1 Hz~ 1 ]. 


(b) The Longitudinal Wave Luminosity 

To obtain explicit forms for the energy fluxes emitted in the form of longitudinal waves, we must 
evaluate the asymptotic Fourier transform at large z (Appendix A), and calculate the product of the indi- 
vidual source functions and their conjugates. We transform the coordinates in equation (3.17) to average 
positions and times of the interacting turbulent eddies and then evaluate the turbulence velocity correla- 
tions by assuming that the fourth-order correlations can be reduced to second-order correlations (cf., 
Batchelor 1953). We calculate the one-dimensional convolutions (Appendix B) over the turbulence 
spectrum by asstiming the spectrum to be a separable product of a frequency-independent energy spec- 
trum and a frequency factor (see section 5a and Appendix B). Finally, after some algebra, we obtain the 
wave luminosity spectrum [ergs s~ l Hz~ l ], 


L t (k,a>j) = 


T„ 


(1 + 5) 


P& k - iX 
©V, 4 k 2 


I S, (k ,m) 1 2 sin 2 (fe - <{>) , 


(4.13) 


where $ is an arbitrary phase (henceforth assumed to be zero). The sin 2 (fcr - $) term appears as a result 
of squeezing the tube at regular intervals of z + nrc/k along the tube; at these points, the efficiency of 
longitudinal wave generation is highest. In addition, I £,(£,©) 1 2 is the product of the individual source 
function and its conjugate, and is given (from 3.17) by 


IS,(*,©)I 2 = j- p 2 / e (fc,©) 


(P?© 2 + ©S) 2 + -£B, 4 c o 2 k 2 V, 2 


_16 

f 


and where the convolution integral /„.<*,©) is defined (Appendix B) by 


(4.14) 


/ e (*,©) = -±r j dk'E 3 (k')E 2 (q)g(kr,( 0 ), 


(4.15) 
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where k' is the wave number of a turbulent eddy, and q = k - k'. As is customary, E 2 (q t ) and E^q) 
denote the one-dimensional turbulence energy spectra evaluated at the two distinct wave numbers q, and 
q. Following Heinze (1975), we note that the one-dimensional (1-D) energy spectrum E^q) is con- 
nected to the three-dimensional (3-D) energy spectrum £ (fc 0 by the relation 


Ei(.q) = dk'£p- 


1 + 


sL 

k' 2 


(4.16) 


where k' is the wave number of all turbulent motions which contribute to the 1-D energy spectrum. In 
addition, the function g(kjt',<&) is expressed by die frequency factor. Explicit forms for both the energy 
spectrum E(k 0 and the frequency factor are given in subsection 5a below. 

The imaginary term in equation (4.13) can be neglected, as it vanishes upon integration over eo. In 
addition, we spatially average L,(zj) over height < Uk^ by performing a spatial integration of equa- 
tion (4.13), and obtain die height-dependent contribution to the monochromatic wave luminosity 
[ergs s~ l Hz' 1 ] for a given flux tube. 


L f> , a (fc,{ 0,z ) a «V0>, # = 

/ e (*.co) sin 2 (fe) , (4.17) 

where P 2 = (^/(oj 2 - Q, 2 ), and a factor of 1/2 was included to take into account the fact that only the out- 
going flux is considered. 

Separating the dimensional factors by using the turbulent velocity u, and the turbulent length scale 
l,, and performing the integration over r and a>, we obtain the total luminosity [ergs r -1 ] due to longi- 
tudenal flux tube waves 

H 

L * N, Jdz Jd<B ,0)) = 

H 

= 4tc N, f dz 


J,A 0 £,<b, (1 + 5) Af, 3 | dm p£p e (a* + dj) J e (k,m) si n?(*z) , (4.18) 


= 2jl4o (1 + 5) 


P«P,P&Pc 

V, 


(Pa tO 2 + C O 2 ) 


2 ' 2 16p> 2 


fP 2 K, 2 


where N, is the number of flux tubes on the stellar surface, which can be related to the magnetic flux 
tube surface filling factor E, ( = p,u, 2 /2) is the kinetic energy of turbulent motions; co, (= «,//,) is the 
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characteristic turbulent frequency; M, ( = u,/V t ) is a coupling Mach number, and the convolution integral 
J c (kM, k, and 55 are dimensionless. Note further that H is the thickness of the turbulent region in the 
convection zone, and that in the derivation of equation (4.5a), we have used the dispersion relation to 
calculate the wave number. In addition, a* and a* are the coefficients of monopole and dipole emis- 
sion, respectively, and are defined by 


and 


Om 





4P, 4 

(Y~ 1) Pi 


(4.19) 


a* ■ ® 2 p, 4 , (4.20) 

where J5 CI is the dimensionless critical frequency (equation 3.4). 

Note that both emission coefficients are calculated by including the term UaJc (which comes from 
the evaluation of the asymptotic Fourier transform, see Appendix A) in the source function (equation 
4.14), and by separating all terms that show dependence on ffi 2 , and on ffi 2 ; the former and latter terms, 
respectively, are the coefficients of monopole and dipole emissions with respect to the wave frequency 
GS. Note that because of this procedure, the coefficients of monopole and dipole emission, given by 
equations (4.19) and (4.20), are redefined in comparison to equation (3.13), and that now only the time- 
dependent term (which derives from non-stationary turbulence) accounts for the coefficient of dipole 
emission; the latter emission can be neglected if stationary or quasi-stationary turbulence is assumed. 
Both emission coefficients are dimensionless, and depend on the turbulent energy spectrum, the physical 
conditions in the region where the waves are generated, and on the critical frequencies for longitudinal 
waves. 

The total wave luminosity for longitudinal tube waves (equation 4.18) shows a dependence on the 
third power of the Mach number (dipole type of emission with respect to Mach number); this distin- 
guishes our results from those obtained by Stein (1968), who considered the generation of acoustic 
waves in a non-magnetic stratified medium, and found dependence of the acoustic wave luminosity on 
the fifth power of the Mach number (quadrupole type of emission). The results presented here are also 
different from those given by Musielak and Rosner (1987a), who obtained monopole type of emission 
for the generation of compressible MHD slow waves in a stratified medium with an embedded uniform 
magnetic field. 

To calculate the energy fluxes given by equations (4.18) - (4.20), we have to specify the number of 
flux tubes on the stellar surface, as well as die values of the magnetic field strength, pressure and den- 
sity at the level in the atmosphere where the integrations are to take place. In addition, we must deter- 
mine the wave frequency domain for longitudinal waves, and must describe the turbulence; in the latter 
case, we must know the shape of the turbulence energy spectrum, the frequency factor, and the 
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characteristic length scale of the turbulent motions. In the next section, we will show how these param- 
eters are restricted by the observational data or by theoretical studies. 
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5. MODEL PARAMETER DEPENDENCE OF ENERGY FLUXES 

In this section, we present preliminary energy fluxes and energy spectra for longitudinal waves, 
and discuss the dependence of our results on the free parameters. We begin with a description of the 
turbulence, discuss the dependence of the energy fluxes on the magnetic field and, finally, present the 
energy spectra for longitudinal tube waves and show how the results depend on the chosen magnetic 
field strength, the turbulent energy spectrum and the frequency factor. 

(a) The Turbulence Energy Spectra and Turbulent Length Scale 

Turbulent motions are characterized by the turbulent velocity and by the turbulent length scale; 
both these parameters are necessary to calculate the turbulent energy spectra, but have a simple form 
only for the very special case of isotropic, homogeneous and incompressible turbulence. In this simple 
case, dimensional analysis shows that the frequency-independent energy spectrum E(Jc) has the Kolmo- 
gorov form 



(* **,) 


(5.1) 


in the inertial range; alternatively, at smaller scales, where viscous effects begin to play a role, the spec- 
trum takes on the exponential form 


E{k) = 


6W k 
k, k, 


| exp[-4*/* t ] , 


where the normalization condition for E(k) is given by 


(5.2) 


[Eqt)dk = | u, 2 , (5.3) 

and where the ma xima of the Kolmogorov and exponential spectra occur near k, = lull,. 

For the frequency factor, either the Gaussian form 


id 12 ku t 


A(*,(o) = 


exp [-{(Orta*) 2 ] 


(5.4) 
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or the exponential foim 


A (k ,co) = exp [-(co/faik) j , 


can be assumed (see Stein 1967), with 


«* 



(5.5) 


(5.6) 


In addition, A(£,co) is normalized so that 


I A(*,to)dco = 1 . (5.7) 

The shape of the turbulence energy spectrum for stratified atmospheres is not known, and neither 
observational data nor theoretical considerations can help us to deduce its general properties. The 
energy spectra presented above certainly do not apply to the stratified turbulent atmosphere, especially to 
the largest turbulent eddies, whose size is comparable to the atmospheric scale height. Nevertheless, we 
shall adopt the above forms for the spectra since they span the range of likely behaviors of the actual 
spectra. 

As an aside, we note that Stein (1967, 1968) and Bohn (1980) similarly adopted various ad hoc 
forms for the energy spectra, using the two functional forms just discussed and the so-called Spiegel tur- 
bulent spectrum, and three different functional forms for the frequency factors; in all these cases, strati- 
fied and magnetic field-free atmospheres were considered, and a turbulence correlation length scale equal 
to the pressure scale height was assumed. In addition, these authors assumed that the turbulent velocity 
of the largest eddies is the same as the velocity of convection motions given by the convection zone 
model In our approach, the turbulent medium which surrounds a tube is stratified and magnetic field- 
free, so that in order to describe the turbulence, we have to make the same assumptions as those just 
mentioned. The difference is, however, that the largest turbulent eddies (with sizes comparable to the 
density or pressure scale height) do not contribute to the generation of longitudinal tube waves; in the 
interaction between the tube and turbulence, only eddies with sizes comparable to the tube diameter are 
important, and these eddies dominate longitudinal wave generatioa We will include this effect in future 
energy flux calculations by computing the diameter of tube d, for a given height (using the magnetic 
flux conservation law) and then, by setting l, = d,, estimate the turbulent velocity of the eddies with size 
d t . Unfortunately, neither the variations of the tube diameter with height in the deep photospheric 
layers, nor even the diameter of an individual "elementary” flux tube can be estimated with any confi- 
dence from the observational data because of current limitations on spatial resolution (cf., Solanki and 
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Stenflo 1985). 


(b) The Magnetic Field Strength 

In the approach presented in this paper, the magnetic field strength is a free parameter, but its max- 
imum value is in fact fairly constrained by observational data, as well as by the horizontal pressure bal- 
ance of the flux tube (see section 2d). From the observational point of view, the magnetic field strength 
within flux tubes can be estimated for the Sun by analysing the statistical properties of the Stokes I and 
V line profiles (Solanki and Stenflo 1985); this method gives the field strengths in typical network 
regions ranging from 1,400 G to 1,700 G. From a more theoretical perspective, horizontal pressure bal- 
ance (based on the thin flux tube approximation) restricts the magnetic field pressure within the tube to 
be less than or equal to the gas pressure outside the tube (equality only obtains in the unlikely case that 
the flux tube is purely "magnetic", i.e., that there is no gas inside the tube). For typical published 
models erf the solar photosphere (for example, Vemazza et al. 1983), the maximum field strength 
obtained from horizontal pressure balance also does not exceed 1,700 G. 

In our approach, the horizontal balance for pressure and magnetic field perturbations (equation 
2.26) restrict p x and B x to be quantities of the same order as the turbulent pressure outside the tube; note 
that the turbulent pressure is a small fraction of the gas pressure outside the tube. If the gas pressure 
inside the tube becomes comparable to the turbulent pressure outside the tube, then perturbed and unper- 
turbed quantities become comparable and the inhomogeneous tube wave equation (3.1) is no longer 
valid. In order to insure the validity of the perturbation scheme, we therefore consider only the cases 
when the unperturbed gas pressure inside the tube is at least twice as large as the turbulent pressure 
defined by equation (2.25); our method of calculation does not allow us to estimate the energy fluxes 
generated in the form of longitudinal waves for magnetic field strength values close to the maximum 
value. 


(c) The Wave Luminosity Spectra 

The wave luminosity spectra presented in this paper are preliminary, and are obtained for log g - 
4.5, for one fixed tube diameter l t =d,= 0.5 H r for one fixed Mach number M = u,IV, = 0.1 and for 
magnetic field strengths varying from 1,000 G to 1,550 G; the latter variations of the magnetic field lead 
to variations of the tube Mach number from 0.15 to 0.08, respectively. These assumptions significantly 
simplify the problem, and allow us to show dependence of the wave luminosity spectra on the chosen 
magnetic field, the turbulent energy spectrum and the frequency factor. Note, however, that we do not 
need to specify the number of flux tubes on the stellar surface (or the filling factor) as we do not calcu- 
late the total wave luminosity. 

The contribution of monopole and dipole emission to the total wave luminosity spectrum is shown 
in Figure 1. The results presented in this figure are obtained for the exponential energy spectrum and 
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the exponential frequency factor, and for values of the magnetic field strength of B„ = 1,500 G and 
1,100 G. In the both cases, monopole emission show a maximum at the same wave frequency 
( co £ 5£l t ) as dipole emission; however, the efficiency of wave generation by dipole emission always 
exceed that of monopole emission; the latter conclusion is valid for both turbulent energy spectra and 
both frequency factors considered here. However, when the magnetic field strength decreases, the effi- 
ciency of monopole emission increases, and becomes comparable to that of dipole emission for B 0 < 
1,200 G; this occurs mainly because P, increases when the field strength decreases. 

The dependence of the frequency-integrated energy flux on the magnetic field strength is given in 
Figure 2, which shows the results obtained for the exponential turbulent energy spectrum and the 
exponential frequency factor. The results strongly depend on the magnetic field strength: the energy flux 
decreases when the magnetic field strength increases. The latter effect leads to a decrease of the gas 
pressure inside the tube, and thus to a decrease of the wave energy flux (as there is not enough gas pres- 
sure to support the waves); this effect also places a limit on the longitudinal wave generation rate. The 
ma ximum of the wave luminosity spectrum occurs well above the cutoff frequency (5 £1,). 

Figure 3 shows the dependence of the wave luminosity energy spectra on the form of the chosen 
turbulent energy spectrum and the frequency factor. As discussed above, we considered Kolmogorov 
and exponential energy spectra (equations 5.1 and 5.2) and Gaussian and exponential frequency factors 
(equations 5.4 and 5.5). The results show that the efficiency of wave generation is higher for the 
exponential energy spectrum, and also that in this case the maximum of the spectrum occurs at higher 
frequencies. 
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6. DISCUSSION AND CONCLUSIONS 


We have derived the source function and the energy fluxes for the generation of longitudinal tube 
waves in a thin, vertical magnetic flux tube embedded in an otherwise magnetic field-free stellar convec- 
tion zone; and we have shown that dipole emission dominates in the generation of longitudinal tube 
waves in such flux tubes; earlier qualitative results given by Stein (1981) and Ulmschneider and Stein 
(1982) suggested dominance of monopole emission. Our results distinguish tube wave generation from 
acoustic and MHD wave generation in a stratified medium since in these latter cases, quadrupole and 
monopole type of emissions are expected to be dominant (Stein 1968, Musielak and Rosner 1987a). 
The efficiency of longitudinal wave generation decreases significantly when the magnetic field strength 
increases, mainly the decreasing gas pressure inside the tube cannot support large energy fluxes for such 
waves. Note also that since the gas pressure within the tube decreases with increasing magnetic field 
strength, our approximations for the generation rate of longitudinal tube waves become unreliable when 
the gas pressure inside the tube is comparable to the turbulent pressure outside the tube. 

When acoustic waves are generated in a stratified medium, monopole and dipole emissions also 
appear, and show a narrow maximum in the energy flux for wave frequencies very close to the acoustic 
cutoff frequency (Stein 1967, Figure 5); this energy flux cannot be carried away efficiently by acoustic 
waves from the generation region because these waves are almost evanescent. The situation is different 
for the generation of longitudinal waves; the maximum of the energy flux generated by dipole emission 
occurs at to = 5fl r , and is somewhat broadened for higher frequency waves (mainly because the source 
function does not depend on the tube cutoff frequency, but instead depends on the Brunt-Vaisala fre- 
quency modified by the tube geometry). The latter frequency is significantly lower than the tube cutoff 
frequency, which is comparable to the acoustic cutoff. Even if the acoustic cutoff frequency co* (as 
well as the tube cutoff frequency £2,) increases in stellar photospheres (for example, to* = 0.024 at the 
bottom of the convection zone and w*. = 0.034 at the temperature minimum; see Vemazza et al. 1976, 
Model C), we may conclude that the maximum of the longitudinal wave energy flux reaches the tem- 
perature minimum region as well as the chromosphere. To estimate however the amount of longitudinal 
wave energy flux that is available for the heating of the lower and upper chromosphere, further calcula- 
tions of longitudinal wave propagation, which take radiative damping into account, are necessary. 

The energy fluxes generated as longitudinal tube waves strongly depend on the strength of the 
magnetic field and the number of flux tubes on the stellar surface; any variations of these parameters 
among late-type stars with the same spectral type will lead to different levels of chromospheric activity 
observed for the same spectral type. Finally, we note that work corresponding to that discussed here 
remains to be carried out for transverse tube waves; this will be the aim of a succeeding paper. 
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APPENDIX A: The Asymptotic Fourier Transform 

Let us consider a one-dimensional inhomogeneous partial differential equation with constant coeffi- 
cients, given in the form 


P 


' jP_ JP_ 

dt 2 ’dz 2 

V 


u(z^) = f(.zj) , 


(Al) 


where P is a polynomial and / (z ,r) is a function which vanishes outside a restricted region. To solve 
this equation, we Fourier-transform u(zj) and f(zj), substitute these into equation (Al), and obtain 


!/(*,©) = 


F(k, co) 

P (-to 2 , -A 2 ) ’ 


(A2) 


where U (k ,co) and F(k, co) are the Fourier transforms of u(zj) and f(zj), respectively. The solution for 
u(zj) can be written as 


u(z,r) = J das J dk exp[-i(cor - kz)) — - ; 


(A3) 


for the special case of a monochromatic source with frequency co 0 , 


we have the result 


F(*,(o) = F{k) 5(oo - co,,) , 




(A4) 


(A5) 


Dropping the subscript o , and since the integration over k can be performed with the help of the residue 
theorem, we find 


u(zj) = exp(-<D0 2jri 


exp (ik H z)F(kn) 

t iGCMo) 

At A 


(A6) 
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where G(jfc,co) = P ( -to 2 , -k 2 ); see also Lighthill (1960). 

Using equation (A6), we may evaluate the asymptotic Fourier transform of /„( 2 .to) given by equa- 
tion (4.11), and obtain 


/„(*,cv) = - ^ (k -iX) ^S*(£,co) txp(ikz) - S, (k ,co) exp(-ifa ) j = 


= 4JI 2 ~ k2 lX l5,(*i,co)l 2 sm 2 (*,z - <p) , 


(A7) 


where 4> is an arbitrary phase. 



APPENDIX B: The Convolution Integral 


In order to calculate the explicit form of the source function, it is necessary to evaluate a fourth- 
order velocity correlation. As this cannot as yet be done from first principles, we follow customary pro- 
cedures (in the present context, see for example Stein 1967) and replace the fourth-order correlation by a 
sum of products of second-order correlations: 


<Ot.it, + UyUy) (it, 'it,' + = 8<it,u x ><it,V> . (Bl) 

where we have used the fact that there is no difference between the x and y directions. We thus have to 
evaluate the convolution integral 

Jxxxx(k,(a) = -“-r j dr j dx <Ux(x y j j)Ux(x y j + rj + x)> x 

( 2 it)* jL 

x <u x Ocyj,t)u x (xy f z + rj + t )> exp[/ (cor - kr)] = 


= J dk’J da>'4> a (k - k\ to - (o') ^(jfc'.co') , 


(B2) 


where (k ,co) is the Fourier transform of the velocity correlations. 

We assume that d>„ (k ,co) may be factored into the frequency-independent one-dimensional energy 
spectrum E 2 (k) and the frequency factor A(i ,co) 


<t>„(*,<o) = -*• E 2 (*) A(*,o>) , (B3) 

where k and u k are the wave number and the velocity perturbation of the k* eddy, and where the factor 
1/4 comes from the different normalization of <&„(*, co) as compared to the energy spectrum E 2 {k) and 
the frequency factor 8. Note that Erfk) is normalized so that 
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and is connected to the three-dimensional energy spectrum (see Heinze 1975). 

Using equation (B3), we may evaluate the convolution integral (B2) and obtain 


/ TO (*,a>) = J dk'Etf - k') g{ky,<o) , (B5) 

where 

g(kjt » = } d (o'A(* Jfc ',(0 - (o') 8 (* .(o'). (B6) 


Because there is no difference between the x and y-directions, we have 


Ar*x(*,co) = / w (*,(o) = = /^(fc.co) = 7 c (i,(o) . (B7) 
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FIGURE CAPTIONS 

Figure 1 : We display the efficiency of monopole and dipole emission for two different values of the 
magnetic field strength, B 0 = 1,500 G (solid lines) and 1,100 G (dashed lines). We assume an 
exponential turbulent energy spectrum and an exponential frequency factor. 

Figure 2: The frequency-integrated wave energy flux for longitudinal tube waves is plotted for different 
values of the magnetic field strength. We assume an exponential turbulent energy spectrum and 
an exponential frequency factor. 

Figure 3: The wave luminosity spectra obtained for a Kolmogorov energy spectrum and an exponential 
frequency factor (KE) are compared to those obtained for an exponential turbulent energy spec- 
trum and an exponential frequency factor (EE), and an exponential energy spectrum and a Gaus- 
sian frequency factor (EG). The magnetic field strength is assumed to be B„ = 1,500 G. 
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